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PREFACE 


This  report  describes  work  done  In  the  summer  of  1983  by  Dr  B.  D. 
Slvazllan,  Department  of  Industrial  and  Systems  Engineering,  the  University  of 
Florida,  Gainesville,  Florida  32611  under  Contract  No.  F08635-83-C-0202  with 
the  Air  Force  Armament  Laboratory  (AFATL),  Armament  Division,  Eglin  Air  Force 
Base,  Florida  32542.  The  program  manager  was  Mr.  Daniel  A.  Mclnnis  (DLYW). 

The  work  was  initiated  under  a  1982  USAF-SCEEE  Summer  Faculty  Research 
Program  sponsored  by  the  Air  Force  Office  of  Scientific  Research  conducted  by 
the  Southeastern  Center  for  Electrical  Engineering  Education  under  Contract 
No.  F49620-82-C-0035. 

This  work  addresses  itself  to  the  problem  of  computing  the  uncertainty 
associated  with  the  probability  of  kill  when  using  the  two-parameter  Carleton 
damage  function  as  specified  in  the  Joint  Munitions  Effectiveness  Manual/Air- 
to-Surface  ( JMEM/AS)  open  end  methods  as  described  in  61  JTCG/ME-3-7  (Revised 
15  May  1980). 

The  two-parameter  Carleton  damage  function  approximates  the  probability 
of  kill  due  to  fragmentation  of  an  exploding  weapon  in  the  absence  of  blast 
effect  and  delivery  error.  Further,  it  assumes  that  at  the  center  of  the 
exploding  weapon,  the  probability  of  kill  is  unity.  It  thus  excludes 
weapon/target  situations  in  which  such  probability  of  kill  results  in  a  number 
less  than  unity. 

The  author  has  benefited  from  helpful  discussions  with  several  people. 
Particular  thanks  are  due  to  Mr  Jerry  Bass,  Mr  Daniel  Mclnnis,  Mr  Charles 
Reynolds,  and  Ms  Katherine  H.  Douglas,  all  from  DLYW  who  have  read  the  report 
and  have  contributed  to  it  through  helpful  comments. 

The  report  is  the  first  of  a  series  dealing  with  the  uncertainty 
associated  with  various  weapon  effectiveness  indices,  and  details 


methodologies  and  techniques  used  in  computing  such  uncertainties  In  the 
presence  of  error  In  the  Input  parameters. 
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SECTION  1 


INTRODUCTION 

The  present  report  discusses  methodologies  for  estimating  the  probability 
of  kill,  P^,  taken  as  a  measure  of  weapons  effectiveness.  In  addition,  the 
report  addresses  itself  to  the  important  question  of  estimating  the 
variability  in  Pk  in  the  presence  of  uncertainty. 

Section  II  elaborates  on  the  techniques  available  to  determine  the  mean 
and  the  variance  of  Pk,  denoted,  respectively,  by  E[Pkl  and  Var[Pk],  when  Pk 
is  expressed  as  a  mathematical  function  of  given  input  parameters,  say  X  and 
Y.  It  is  assumed  here  that  none  of  X  or  Y  are  known  precisely  but  are  subject 
to  estimation  error.  Such  error  may  arise  if,  for  example,  X  and  Y  are 
measured  subjectively,  or  X  and  Y  are  obtained  through  some  type  of 
inferential  estimation  procedure  such  as  the  use  of  multiple  regression 
scheme. 

Section  III  applies  the  methodology  to  the  particular  situation  of 
fragment  sensitive  targets  in  the  absence  of  blast  and  aiming  error.  The 
probability  of  kill  in  such  an  instance  can  be  approximated  by  the  Carleton 
damage  function 

Pkf  -  exp[-(4+4)l  • 
x  y 

Here  it  is  assumed  that  the  weapon  explodes  at  (0,0)  and  the  target  is  located 
at  (x,y).  R  and  R  are  two  parameters  identifying  the  weapon's  radii.  The 
general  approach  assumes  that  estimates  of  Rx  and  R^,  as  well  as  the  errors  in 
such  estimates,  are  available,  and  the  problem  addresses  itself  to  estimating 
?uf  and  its  variance. 


Section  IV  goes  one  step  beyond  and  looks  at  the  problem  of  actually 
estimating  Rx  and  Ry  as  well  as  their  covariance  matrix.  It  is  assumed  that 
experimental  results  obtained  from  fragmentation  field  data  of  an  exploding 
weapon  are  used  to  estimate  the  values  of  Rx  and  Ry  in  the  Carleton  damage 
function.  Two  methods  are  explored,  both  based  on  the  theory  of  linear 
multiple  regression. 

Some  concluding  statements  are  made  in  Section  V. 

At  this  point  two  remarks  are  in  order.  The  first  remark  concerns  the 
Carleton  damage  function.  In  the  cases  considered  in  this  report ,  it  is 
implicitly  assumed  that,  at  the  point  of  weapon  explosion,  the  probability  of 
kill  is  unity.  Further,  one  assumes  that  the  equiprobabil ity  contour  lines  in 
the  (x,y)  plane  are  ellipses  with  axes  coinciding  with  the  x  and  y  axes. 

These  conditions  are  fairly  reasonably  satisfied  when  the  impact  angle  of  the 
weapon  is  close  to  90  degrees  and  the  x  and  y  axes  are  taken  to  be  the 
directions  of  the  weapon  range  and  deflection,  respectively.  Thus,  one  may 
use  the  two-parameter  Carleton  damage  function  as  a  model  for  this 
situation.  When  the  impact  angle  is  smaller  than  90  degrees,  the  equiproba¬ 
bil  ity  contour  lines  do  not,  in  general,  follow  elliptic  patterns  and  are  not 
in  general  symmetrical  about  the  deflection  axis.  However,  because  of  its 
simplicity,  the  Carleton  damage  function  is  still  used  in  practice. 

When  the  probability  of  kill  at  the  point  of  impact  is  less  than  unity, 
the  three-parameter  Carleton  damage  function  given  by 

pkf  *  °o 

x  y 

is  a  suggested  model  (see  e.g.,  [6]).  Notice  that  in  this  case,  since 
0  <  <_  1,  one  must  necessarily  have  0  <  D0  £  1.  In  what  follows,  only  the 

two-parameter  Carleton  damage  function  is  considered  (D0=l). 


The  second  remark  concerns  the  fragmentation  field  data  for  computing 
It  will  be  assumed  that  such  data  are  available  and  that  they  have  been 
processed  to  exhibit  in  matrix  form  the  value  of  at  a  target  point  (x,y), 

when  the  weapon  impacts  at  (0,0)  with  a  given  angle  and  at  a  given  velocity. 


SECTION  II 


METHODOLOGY 


1.  Background 

In  many  practical  situations,  a  measure  of  the  effectiveness  of  an 
exploding  weapon  upon  a  target  is  usually  expressed  by  the  probability  of  kill 
Pk.  The  computation  of  this  Pk  value  is  most  often  undertaken  by  using  a 
mathematical  expression  or  model  which  is  a  function  of  the  target  position 
relative  to  the  center  of  explosion,  as  well  as  a  function  of  one  or  more 
parameters.  These  parameters  are,  in  general,  established  either  subjectively 
or  experimentally  or  both,  and  they  account  for  a  number  of  factors  which  may 
be  present  at  the  time  of  the  explosion,  such  as  weapon  characteristics, 
impact  angle,  impact  velocity,  nature  of  target,  etc.  An  example  of  a  Pk 
function  is  the  Carleton  damage  function  which  was  previously  introduced. 

Thus,  given  a  set  of  conditions,  the  formula  allows  one  to  compute  a  Pk  value 
for  any  location  of  the  target  relative  to  the  weapon.  However,  one  must 
realize  that  this  formula,  like  any  other,  is  only  an  approximation  whose 
usefulness  is  dictated  by  how  accurately  it  represents  reality.  This  is  due 
to  the  fact  that  (a)  the  mathematical  model  is  not  an  exact  replica  of  the 
actual  situation  and  (b)  the  parameters  in  the  model  are  estimates  (random 
variables)  rather  than  exact  quantities. 

A  natural  question  that  may  be  raised  in  this  context  is  how  good  is  the 
Pk  value  computed  from  the  mathematical  model;  or,  precisely  what  is  the  error 
of  the  Pk  value?  To  answer  this  question,  one  should  note  that  Pk  being  a 
function  of  the  estimates  is  itself  a  random  variable.  Theoretically 
speaking,  given  the  joint  distribution  of  the  parameters,  the  distribution  of 
Pk  could  be  obtained,  thus  providing  a  means  for  calculating,  for  any 
confidence  level,  interval  estimates  for  Pk«  However,  in  practice  this  is  not 


>1 
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possible  due  to  several  reasons.  First,  from  an  economic  point  of  view, 
inadequate  data  are  available  to  describe  fully  the  joint  distribution  of  the 
parameters.  Second,  Pk  is,  in  general,  a  complex  function  of  the  parameters, 
and  thus,  characterizing  its  distribution  becomes  very  difficult,  if  not 
impossible.  Finally,  the  Intent  is  not  to  set  up  accurate  confidence 
intervals  on  Pk,  but  rather  to  provide  the  decision  maker  with  adequate 
information  on  how  Pk  behaves  statistically  in  the  presence  of  estimation 
error  in  the  input  parameters.  Thus,  very  often  the  parameters  will  often  be 
characterized  by  their  first  two  moments.  Similarly,  the  establishment  of  the 
first  two  moments  of  Pk  will  be  adequate  for  our  purpose. 

For  argument  sake,  suppose  that  the  probability  of  kill,  Pk,  is  a 
function  of  the  two  input  parameters  X  and  Y  so  that 

Pk  =  Pk  (X,Y)  .  (1) 

Suppose  that  each  of  the  parameters  X  and  Y  are  estimates  subject  to  error. 

It  is  required  to  determine  E[Pk]  and  Var[Pk], 

If  the  joint  distribution  of  X  and  Y  is  known  and  is,  say,  fx  y  (x,y), 

then 

E[Pk]  =  /  /  PR  (x,y)  fx>Y  (x,y)  dx  dy 

and  E[Pk]  =  JX^J  P^  (x,y)  f^y  (x,y)  dx  dy  . 

From  these  two  expressions  one  can  obtain 

Var[Pk]  =  E[Pk]  -  EZ[Pk]. 

It  should  be  noted  that  although  it  is  impossible  to  completely 
characterize  the  function  fx>y  (x,y),  nevertheless,  there  exists  methodologies 
to  arrive  at  approximate  estimates  of  E[Pk]  and  Var[Pk].  We  shall  digress  on 
two  such  methods  which  are: 


a>  the  subjective  estimation  procedure 
b.  the  Taylor's  series  estimation  procedure. 


Each  of  these  methods  assumes  that  certain  statistical  information  is 
available  on  all  input  parameters  which  enter  in  the  computation  of  P^. 

2.  The  Subjective  Estimation  Procedure 

In  the  subjective  estimation  procedure,  it  is  assumed  that  the 
uncertainty  level  of  each  input  parameter  is  provided  as  subjective 
information.  A  lower  and  upper  bound  value  for  each  parameter  is  obtained. 

The  value  of  a  particular  parameter  is  assumed  to  take  equally  likely  values 
between  its  two  extreme  points.  This  is  equivalent  to  assuming  that  each 
parameter  is  a  random  variable  uniformly  distributed  over  its  range  of 
values.  Further,  the  parameters  are  assumed  to  be  mutually  independent  random 
variables.  With  this  statistical  information,  the  evaluation  of  E[Pk]  and 
E[P^]  are  reduced  to  the  computation  of  a  set  of  definite  integrals.  Thus, 
referring  to  (1),  if  X  is  assumed  to  be  uniformly  distributed  in  the  interval 
[xj.x^],  while  Y  is  assumed  to  be  uniformly  distributed  in  the  interval 
Cylty2],  then,  if  X  and  Y  are  mutually  independent,  it  follows  that 


E[Pk]  - 

E[P2]  = 


1 

(x2-x1)(y2-y1) 

1 

(x2-x1)(y2-y1) 


y9  x? 
I  2f  2 
7  xi 


y2  *2 


I  7 


7  xi 


(x,y)  dx  dy 


( x  ,y)  dx  dy  . 


Note  here  that  under  these  assumptions,  the  mean  and  variance  of  X  and  Y  are 
given,  respectively,  by 
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;  Var[X] 


ECX]  = 

ECY]  • 


Xl+X2 


V^2 

2 


(x2"xl^2 

=  — T2 — 


Var[Y] 


[Mill 
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Although  the  method  often  provides  expressions  for  E[Pk]  and  E[Pk]  which 
can  be  analytically  manipulated  to  arrive  at  closed  form  expressions,  never¬ 
theless,  it  has  certain  inherent  disadvantages  which  should  be  stated  at  this 
stage.  First,  there  is  no  guarantee  that  closed  form  expressions  can  be 
obtained  for  E[Pk]  and  E[Pk].  If  such  expressions  are  derivable,  they  are 
usually  fairly  complex  in  form.  Second,  the  method  does  not  make  any 
allowance  for  incorporating  a  dependency  factor  between  X  and  Y,  if  such 
dependency  is  known  to  exist.  A  third  disadvantage  of  the  method  lies  in  the 
fact  that  it  is  not  possible  to  segregate  the  contribution  of  the  variance  of 
each  parameter  component  to  the  variance  of  the  probability  of  kill  Pk. 
Finally,  the  numerical  computation  of  Var[Pkl  requires  that  the  computed 
values  of  E[Pk]  and  E[P^]  be  carried  to  several  significant  digits. 

3.  The  Taylor's  Series  Estimation  Procedure 

In  the  Taylor's  series  estimation  procedure,  the  assumption  is  made  that 
the  statistical  moments  of  the  input  parameters  are  known  and  that  Pk  (X,Y)  is 
a  differentiable  function  of  X  and  Y.  The  expression  for  the  probability  of 
kill  Pk  is  expanded  as  a  Taylor's  series  about  the  expected  value  of  the  input 
parameters.  Only  first  order  terms  are  assumed  significant  in  the  derivation 
that  follows.  However,  in  general,  one  could  use  a  procedure  parallel  to  that 
outlined  here  if  higher  order  terms  beyond  the  first  are  to  be  included. 

Let  J  =  E[X]  and  7  =  E[Y],  then 

_  3Pk  I  -  3Pk  I 

Pk  (X.Y)  “  Pk  (X,Y)  +  (X-X)  -r—-  I _ +  (Y-Y)  I _ .  (2) 


Taking  expectations  on  both  sides  of  (2)  yields 


E[Pk  (X,Y)]  *  Pk  (X,7)  . 


Thus,  in  the  Taylor's  series  estimation  procedure,  the  probability  of  kill 
computed  at  the  mean  values  of  the  parameters  is  taken  as  an  estimate  of  the 
mean  of  the  probability  of  kill.  This  estimate  is  an  approximation  which  is 
adequate  for  most  practical  problems.  Further  accuracy  may  be  obtained  by  the 
inclusion  of  higher  order  terms  in  the  Taylor's  series  expansion.  For 
example,  the  addition  of  the  second  order  term  to  (2)  yields 


1 

2! 


2  ^k 

{(X-T)2_i 

ax^ 


+  2(X-7)(y-7) 


32Pl 


7,7 


3X3Y 


32pk 

♦  (V.T) 

X.T  3Y^ 


X,Y 


which  upon  taking  expectations  gives 


\  Var[X] 


Cov[X,Y] 


82Pl 


3X3Y 


X»Y 


+  |  Var[Y] 


This  last  expression  would  be  added  to  Pk{7,7)  to  improve  the  accuracy  of  the 
value  of  E[Pk(X,Y)]. 

To  obtain  an  estimate  of  the  variance  of  Pk ,  write  (2)  as 


3Pl 


Pk(X,Y)  -  Pk(X,7)  .  (X-7) 


7,7 


(Y-7) 


3Pk 

IT 


X,  Y 


Squaring  both  sides  yields 


3P  c 

[pk ( X » Y )  -  Pk(X,Y)]2  =  (X-X)2  (-^) 


7,7 


2  3Pk 

♦  (Y-T)z  (^f) 


7,7 


+  2(X-7)(Y-7)  (T5r)(-g?i) 


7,7 


(3) 
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Taking  expectations  on  both  sides  of  (3)  and  remembering  that 
Pk(X,Y)  -  E[Pk(X,Y)]  yields 

3P  ^  jp  ^ 

Var[Pk(X,Y)]  =  Var[X]  (-^)  +  Var[Y](-^)  __ 

X»Y  X»Y 

♦  2  CO.CX.Y]  C^rK-jr)  J  Y  ■  W 

Formula  (4)  provides  a  mean  for  computing  the  uncertainty  in  the  value  of 
Pk  (x,Y)  as  a  function  of  the  uncertainty  in  the  values  of  the  input 
parameters  X  and  Y,  namely  Var[X],  Var[Y]  and  Cov[X,Y].  Of  course,  in  case 
when  X  and  Y  are  independently  distributed,  the  covariance  term  vanishes  and 
expression  (4)  involves  only  the  variance  components. 

The  advantages  of  this  method  are  three  fold.  First,  it  is  possible  to 
obtain  fairly  simple  expressions  for  E[Pk]  and  Var[Pk],  Second,  in  case  when 
the  input  parameters  are  correlated,  the  uncertainty  in  Pk  can  reflect  the 
extent  of  this  correlation  through  a  covariance  term.  Finally,  the 
contribution  to  the  Pk  variance  of  each  variance  component  can  be  identified 
and  segregated.  With  an  objective  towards  reducing  the  Pk  variance,  the 
methodology  allows  one  to  breakdown  the  Pk  variance  into  its  components  and  to 
identify  those  input  parameters  with  the  largest  variance  contribution.  It 
should  be  noted  that  the  method  fails  in  cases  where  Pk  is  not  a 
differentiable  function  of  the  input  parameters. 

As  was  mentioned  earlier,  when  using  the  Taylor's  series  estimation 
procedure,  a  decision  has  to  be  made  on  how  many  terms  are  to  be  retained  in 
the  expansion.  As  a  first  approximation,  only  first  order  terms  are  usually 
retained.  Improved  accuracy  in  both  E[Pk]  and  Var[Pk]  could  be  obtained  by 
the  inclusion  of  second  and  higher  order  terms.  This  however,  would  generate 


cumbersome  mathematical  expressions  and  would  require  a  knowledge  of  the 
higher  moments  of  X  and  Y.  In  what  follows,  only  first  order  terms  are 
considered  in  order  to  maintain  the  simplicity  of  the  expressions  derived.  An 
extensive  discussion  of  the  Taylor's  series  estimation  procedure  Is  Included 
in  [2]. 

4.  Remarks 

a.  Experience  has  shown  that  in  cases  where  X  and  Y  are  independently 
distributed,  the  subjective  estimation  procedure  and  the  Taylor's  series 
estimation  procedure  provide  results  for  E[Pk3  and  Var[Pk]  which  are  in  close 
numerical  agreement.  This  will  be  verified  later  on  in  this  report  when  the 
methods  are  applied  to  fragment  sensitive  targets  in  the  presence  of  no  aiming 
error. 

b.  In  establishing  confidence  Intervals  for  Pk,  a  two  standard 
deviation  («=2)  two-sided  confidence  interval  is  selected.  In  such  a  case, 
Chebyshev's  inequality  guarantees  at  least  a  75  percent  confidence  interval 
since  for  <*=2  we  have 

P{E[Pk]  -  avVa'r[PJT  <  Pk  <  E[Pk]  +  a/Var[PkJ}  >  1  -  -±  =  .75  . 


SECTION  III 

APPLICATION  OF  THE  ESTIMATION  PROCEDURES 


The  Carleton  Damage  Function 


For  fragment  sensitive  targets,  in  the  absence  of  aiming  error,  the 
Carleton  damage  function  may  be  used  as  a  model  to  measure  the  value  of  the 
probability  of  kill  due  to  fragmentation,  Pkf,  at  an  arbitrary  point  (x,y)  in 
the  plane  given  that  the  weapon  bursts  at  the  origin  (0,0).  This  damage 
function  has  the  form 


2  2 

pkf  =  exP  ('(“7  +  ^r)]  (5) 

Rx  Ry 

where  Rx  and  Ry  are  two  parameters  known  as  weapon  radii.  Rx  is  the  weapon 
radius  in  the  direction  of  range  which  is  the  direction  of  the  main  weapon 
axis  (the  x-axis).  Ry  is  the  weapon  radius  in  the  direction  of  deflection 
(y-axis)  perpendicular  to  the  x-axis.  In  general,  Ry  >  Rx. 

It  is  assumed  now  that  the  parameters  Rx  and  Ry  are  not  known  exactly  and 
that  they  are  subject  to  estimation  errors.  Under  these  conditions,  it  is 
required  to  estimate  the  error  in  Pkf.  To  do  this,  one  needs  to  evaluate  both 
E[Pkf]  and  Var[Pkf]  in  order  to  specify  an  interval  estimation  for  Pkf  for  any 
(x,y).  The  estimates  of  Rx  and  Ry  are  assumed  to  be  given  either  in  the  form 
of  a  minimum  value  and  a  maximum  value  through  a  subjective  estimation 
procedure,  or  in  the  form  of  a  vector  mean  and  a  covariance  matrix  through  an 
inferential  estimation  procedure.  Given  the  estimates  of  Rx  and  Ry,  the 
problem  is  to  determine  E[Pkf]  and  Var[Pkf].  We  shall  consider  next  the 
application  of  each  of  the  two  methods  discussed  in  Section  II  to  this 
particular  problem. 


IJ  *  ij»  II  ■  ■  « 


2.  The  Use  of  the  Subjective  Estimation  Procedure 

Here  It  Is  assumed  that  R  and  R  are  uniformly  and  independently 

x  y 

distributed  over  the  respective  ranges  Rx  <  Rx  <  Rx^  and  Ry,  '  Ry  <  V 

The  quantities  Ry  ,  R„  .  Ru  ,  and  Ru  are  supposed  to  be  known.  They  may  be 
X1  x2  yl  “2 

determined,  for  example,  through  a  subjective  procedure  in  which  individuals 
are  requested  to  provide  a  lower  and  upper  bound  on  the  values  of  Rx  and  R^ 
based  on  their  judgement  and  their  experience.  The  main  objective  is  to 
determine  E[Pkf]  and  VarfP^]  . 

a.  Estimation  of  E[P^] 


Using  (5),  the  expectation  of  is  given  by: 


To  proceed  further,  one  needs  to  evaluate  integrals  of  the  form 
B  2 

J(A,B)  =  /  exp(-  —•)  du  . 

A  u* 

Now  making  the  change  in  variable  v  =  ^  in  this  integral  expression  yields 
l 

A  .2  2  . 

J(A,B)  =  /  e"K  v  •  -4  dv  . 

I  v 

B 

Integrating  by  parts  yields 


12 


J(A,B)  =  -ie-‘v  (-!)(- 2k^v)e-KV  dv 

IT 


=  B  exp(-  -  A  exp(-  -  2k2  /  e"k  v  dv  . 

I 


w  1 

In  the  integral  expression  let  kv  =  —  ;  then  dv  =  -  dw  and 


ft 


kv? 


2  2  k  ft  *r 

J( A,B)  =  B  exp  (-  — *■)  -  A  exp(-  — *■)  -  k/2  J  e  2  dw 
B^  A^  k/2 


kv?  w* 


2  2 

=  B  exp (-  JL)  -  A  exp (-  k  )  -  k/?  [/  n  e  fc  dw 

a^  L  n 


B 


0 


kv?  w£ 

J^e'^dw 

0 


kv?  * 


2  2 
B  exp  (-  — y]  -  A  exp(-  — *■]  -  2k /it  [  ■i  f  6  dw 

D  ft  ^  M  A 


B 


0 


kv?  vr 

.  _L  /  B  e  "  2  dw  ] 
fti  0 


Let 


•(*)  -  “ :  /  e 

V?TT  0 


dw 


Thus,  *(z)  is  the  area  under  the  standardized  normal  curve  to  the  right  of  z=0. 
These  values  are  readily  available,  and  can  be  obtained  from  existing  tables 
(see  e.g.,  [1]).  The  expression  for  J(A,B)  can  be  written  as: 


(7) 


The  expression  for  E[Pkf]  given  in  (6)  can  be  written  using  (7)  as 


b.  Estimation  of  Var[Pkf] 

The  following  expression  results  after  squaring  both  sides  of  (5) 

P^f  =  exp[-(~- + -^-J]  .  (9) 

Rx  Ry 

It  is  thus  evident  that  E[Pkf]  can  be  obtained  directly  from  E[Pkf]  given  in 
(8)  by  replacing  x  and  y,  respectively,  by  x/2  and  y/2"  .  This  substitution 
yields 


The  expression  for  Var[Pkf]  can  be  obtained  by  using  (8)  and  (10)  in  the  formul 
Var[Pkf]  =  E[Pkfl  -  (E[Pkf]}2  .  (11) 


c.  Example 

It  is  required  to  determine  an  estimate  of  the  probability  of  kill  Pkf 
for  a  fragment  sensitive  target  located  at  x=50  ft  and  y=100  ft  when  a  weapon 
detonates  at  (0,0).  The  blast  effect  is  neglected,  and  it  is  assumed  that  the 
Carleton  damage  function  given  by  expression  (5)  can  reasonably  be  used  to 
compute  Pkf.  A  subjective  evaluation  of  the  weapon  radii  identifies  the  value 
of  R  to  be  at  least  80  ft  but  no  more  than  90  ft.  Similarly,  the  value  of  R 

A  J 

is  judged  to  be  at  least  160  ft  but  no  more  than  180  ft. 

To  proceed  with  the  calculations,  it  is  assumed  that  the  parameters  Rx 

and  Ry  are  independently  and  uniformly  distributed  over  the  respective  ranges 

80  <  Rx  <  90  and  160  <  Ry  <  180.  Further,  the  requirement  is  to  calculate  not 

only  E[Pkf],  but  also  Var[Pkf]  in  order  to  place  confidence  bounds  on  Pkf. 

Note  that  R„  =  80  ft,  RY  =  90  ft,  Rv  =  160  ft,  and  Rv  =  180  ft. 

X1  x2  yl  y2 

The  expression  for  E[Pkf]  is  given  by  (8).  Substituting  for  the 
numerical  values  of  x,y,  Rx^,  Rx^,  Ry^  and  Ry^  yields 

2  2 

E[pkf]  =  rreo-Tffoy(9ir-w '  f90  exp^"  ©  1  ■ 80  exp^'  ©  J 

-  (2H50K7  [*(^)  -  {180  exp[-  ] 

-  160  exp[ -  (|||)2]  -  (2) (100)  A 

=  ( 20)(To)~  {<*.099,930,47  -  54.130,707,69 

-  100  /7  [*( .883,883,5)  -  *( . 785 ,674 ,2) ]} 

•  (132.199,860,9  -  108.261,415,4  -  200  /?  [*(.883,833,5) 


-  *(.785,674,2)1)  . 


(12) 


The  values  of  the  function  #(  •)  are  obtained  from  the  tabulated  values  in  [1] 
using  linear  interpolation.  As  an  example,  to  compute  #(.883,883,5),  the 
following  two  values  are  read  from  the  table: 

*(.90)  =  .315,939,88 

*(.88)  =  .310,570,35  . 

It  thus  follows  that 


*(.883,883,5)  =  *(.88)  *  E ^ 


=  .310,570,35  + 


.005, 369, 53)(. 003,883, 5) 


=  .311,612,98  . 

One  obtains  similarly 


=  .283,961,44. 

Substituting  for  these  numerical  values  in  (12)  yields 
E[Pkf]  =  .499,582,48. 


In  a  similar  fashion,  using  (10)  one  obtains  for  E[Pkf] 


Etpkfl  =  (180-160) (90-80)  (90  exp t"2  (90 )  1  "  80  exp['2  ^80^  1 

-  ( 2) ( 50)  m 

(180  exp  [-2(ygjj)  -  160  exp[-2(-^)  ] 

-  (2) (100)  /2iir  [»(-(2-)jg§°-)-)  -  *(ii4§ri)]} 

=  (48.546,675,65  -  36.626,668,94 

-  100  /U  [*(1.25)  -  *(1.111,111,1)]} 

{97.093,351,30  -  73.253,337,88} 

-  200  /U  [*(1.25)  -  *(1.111,111,1)]}  . 

From  the  tables  one  obtains  after  using  linear  interpolation 
*(1.25)  =  .394,338,81 
and  *(1.111,111,1)  =  .366,727,93  . 

Substituting  these  numerical  values  in  the  expression  for  E[P2f]  yields 
E[P2f]  =  .249,898,556,3  . 


Using  (11),  the  following  numerical  value  for  VarfP^]  is  obtained 


Var[Pkf]  =  E[Pkf]  -  {E[Pkf]}^ 


=  .249,898,556,3  -  ( .499, 582, 48)2 
=  .000,315,902  . 

The  standard  deviation  of  Pk^  is: 

op  =  /  Var[P  ] 

Kkf  KT 

=  /  .000,315,90?  =  .0178  . 

Thus,  using  Chebyshev's  inequality,  with  at  least  a  75  percent  confidence 
interval ,  we  have 

Pkf  •  E[Pkf]  *  2o 

kf 

=  .4996  _+  .0356  . 

3.  The  Use  of  the  Taylor's  Series  Estimation  Procedure 

Referring  once  more  to  the  Carleton  damage  function  given  by  expression 
(5),  the  parameters  Rx  and  Ry  are  assumed  to  be  subject  to  estimation  error 
and  thus  cannot  be  determined  accurately.  These  parameters  are  to  be 
specified  at  least  in  terms  of  their  first  two  moments,  that  is,  in  terms  of 
their  statistical  means 

Rx  =  E[Rx]  and  Ry  =  E[Ry]  (1 

and  in  terms  of  a  covariance  matrix  defined  by 
/  Var[Rx]  Cov[Rx,Ry]\ 


Cov[R  ,R  ]  Var[R  ] 


(1 


In  case  Ry  and  R  are  subjectively  estimated  in  terms  of  their  ranges,  say 
x  y 

R„  <  R„  <RV  and  Ru  <  R„  <  Ru  ,  and  assuming  that  R„  and  R„  are 

xi  x  x2  yj  y  y2  x  y 

independently  and  uniformly  distributed  over  such  ranges,  the  statistical 
means  and  the  covariance  matrix  are  given,  respectively,  by: 


R  = 
x 


R  +R 
X1  X2 


and 


R  +R 
yl  y2 


l 


0 


(k  -k  y 


12 


(15) 


(15) 


On  the  other  hand,  it  is  possible  that  the  parameters  Rx  and  R  are 

estimated  using  a  specific  inferential  procedure,  such  as  a  linear  multiple 

regression  analysis.  In  such  a  case,  one  may,  for  example,  use  experimental 

data  to  fit  the  Carleton  damage  function,  and  as  a  result  of  such  an  analysis, 

derive  estimates  for  the  first  two  moments  of  Rx  and  Ry  (see  Section  IV). 

We  next  show  how  one  can  use  the  estimates  of  R  and  R  as  expressed  in 

x  y 

(13)  and  (14)  to  derive  expressions  for  E[Pkf]  and  Var[Pkf]. 
a.  Estimation  of  E[Pkf] 

As  a  first  approximation  one  can  write  (see  Section  II). 

*  pkf(«x.  v  • 

Using  expression  (5)  for  the  Carleton  damage  function  function  one  obtains 

2  2 

ECPkf]  -  exp[-(z?  ♦  h)\  •  (17) 

R.  R . 


19 


b.  Estimation  of  Var[P^] 


The  procedure  outlined  in  Section  II  is  used  to  arrive  at  an  expression 

for  Var[Pkf].  Expanding  the  expression  for  P^  =  Pkf  (Rx,Ry)  as  given  in  (5) 

as  a  Taylor’s  series  about  the  point  (R  .RJ  and  retaining  only  first  order 

x  y 

terms  results  in: 


pkf  <W  =  pkf  <R«*y  * (R*  ‘  V  tar1),  . 

X  K  «K 
x*  y 


(Ry  -  Rv)  (aF1) _ • 

y  y  dKv  %  ^ 

J  x*  y 


Transposing  and  squaring  both  sides  results  in 


tpkf  tRx’Ry)  -  Rkf  (Rx-R,)]2  ■  <VR/  tasf),  , 


-  9  3Pkf 

IT  (  Ki'i 


x  ’V 


x  Rx,Ry 


ap.,,  aPi,f  aP|,f 

+  <Ry-Ry)  (.sf),  ,  *  «Rx-Rx)(R,-R,)  l (air-) (557) •  <> 

y  Vy  X  y  *x,  y 

Taking  expectations  on  both  sides  of  (18),  one  obtains  as  a  first  approxima- 


Var[Pkf]  =  Var[Rx]  (3/-)  _  *  Var[R  ]  (p-) 


'x  R  ,R 
x  y 


3Pkf 

J  (sr1).  , 

1  V*, 


3  P  3P 

+  2  Cov[Rx,Ry]  [(3^  )(gp  )]_  —  • 
x  y  Kx,Ky 


Now  from  expression  (5) 


ffl 


y  Ry  Rx  Ry 


Substituting  (20)  and  (21)  in  (19)  yields 


Var[Pkf]  =  {exp[-2(^  +  ^)}\  Var[Rx) 

R  R  R 

x  x  y 

+  (exp[-2(^-  +  ^)]}  Var[R  1 

y  X  y 

,  (exp[-2(4  +  4)]}  Cov[R  R  ]  . 

n n J  nc  A  J 


R  R 
x  y 


R  R 

x  y 


C.  Examc 


Consider  the  data  given  in  the  example  of  Section  III.2.C.  where 
80  <  Rx  <  90,  160  <  Ry  <180,  x  =  50,  y  =  100.  Using  the  expression  (15)  and 
(16)  yields 

».  *R, 

px  «  .  80*%  .  85  ft 

^  ,  160*180  ,  ,;0  ft 


_  W 

Var[Rx]  -  12 


)-80)‘  _  100  Ct,Z 
12  "  12 


{Ry  )  ,  x2 
Var[Ry]  .  JlJ. L_  .  iMjieol!  .  m  ft2# 


By  assumption,  Cov[Rx,Ry]  =  0. 

Using  these  numerical  values  in  (7)  and  in  (22),  we  obtain 


■  «p-!(§)2  *  ($}H 


»  .500 


Var[Pkf]  =  HUM  (eXp[-2(  * 

(85>6  (86)2 


*  liliiMl  ,,xp[.2( 

( 1 70) 6 


(50)2  , 
(85)2 


=  .000,138,4  +  .000,138,4 
=  .000,276,8  . 


The  standard  deviation  of  P^  is 

oD  =  /'"TTT0'fT,'27ff;8  =  .0166  . 

Hkf 

Using  Chebyshev's  inequality  with  at  least  a  75  percent  confidence  interval  we 
have 


pkf  =  EtPkf^  -  2ap  , 
kf 

=  .500  +_  .033. 

This  last  result  is  fairly  close  to  the  numerical  result  in  the  example 
of  Section  III.2.C.  The  computational  simplicity  of  this  last  method  as 
compared  to  the  subjective  estimation  procedure  should  be  pointed  out.  In 
addition,  one  notes  that  the  variance  contribution  of  the  input  parameters  Rx 
and  Ry,  respectively,  to  the  total  Pkf  variance  are  equal.  Of  course,  this  is 
coincidental  and  is  due  to  the  particular  numerical  data  used.  In  general, 
the  variance  contributions  will  be  different.  The  methodology  provides  a 
means  for  segregating  the  variance  component  of  each  of  the  input  parameters 
entering  in  the  computation  of  and  from  that  point  of  view  is  definitely 


more  advantageous  than  the  subjective  estimation  procedure. 


SECTION  IV 


INFERENTIAL  ESTIMATION  PROCEDURES 


1.  Background 


The  next  logical  question  that  arises  consists  in  the  way  the  statistical 
characteristics  of  the  input  parameters  Ry  and  Ru  are  established.  It  was 

a  y 

seen  that  one  of  the  methods  is  the  so-called  subjective  method.  The  other 
method  is  the  inferential  one  in  which  experimental  data  are  used  to  estimate 
the  values  of  the  parameters  Rx  and  when  is  expressed  as  a  mathematical 
function  of  these  parameters.  Two  regression  schemes  will  be  discussed 
through  which  one  could  compute  such  values  as  E[Ry],  E[R  ],  Var[R  ],  Var[R  1 
and  Cov[Rx,Ry],  all  of  which  are  used  in  the  Taylor's  series  estimation 
procedure  to  arrive  at  values  of  E[Pkf]  and  Var[Pk^]. 


2.  A  Linear  Multiple  Regression  Scheme:  Method  1 

The  starting  point  is  the  Carleton  damage  function  which  was  defined  by 
expression  (5)  and  which  is  repeated  here: 


2  2 

Pkf  =  exp[-(^-  +  ^)]  .  (23) 

R  R 
x  y 

It  is  assumed  here  that,  through  the  analysis  of  the  fragmentation  data,  a 
value  of  Pkf  is  specified  at  a  given  location  (x,y)  when  the  weapon  bursts  at 
(0,0).  The  main  objective  is  to  estimate  the  values  of  Rx  and  Ry. 

Suppose  that  there  are  n  data  points.  For  the  ith  data  point 
( i*l ,2,. . . ,n)  Pk^  represents  the  probability  of  kill  due  to  fragmentation  at 
the  point  (x^,y^).  In  what  follows  it  shall  always  be  assumed  that 
0  <  Pkf.  <_  1  for  all  i.  In  case  Pk^  =  0  such  data  point  is  discarded. 


Taking  logarithms  on  both  sides  of  (23)  yields 
-In  Pkf  =  ln(p~)  =  ^  +  *2  * 


kf  RJ  Rj 


In  (24)  let 


z  =  ln(~),  u  =  x2,  v  =  y2 


a  •  -4  and  b  =  . 

R^  R^ 

x  y 


Substituting  these  values  in  (23)  one  obtains 


z  *  au  +  bv. 


For  the  ith  data  point,  we  shall  assume  that 


ei  =  zi  -  au..  -  bvi  ( 1«1.2 . n) 


are  independently  distributed,  and  further  that 


E[ei]  =  0 


Var[e^]  =  a  . 


The  parameters  a,  b,  and  o  are  to  be  estimated. 


At  this  stage  one  may  be  tempted  to  assume  that  z  ■  ln(= — )  Is  normally 

2  *kf 

distributed  with  mean  (au+bv)  and  variance  c  ,  and  attempt  to  use  the  maximum 

2 

likelihood  technique  to  arrive  at  estimates  of  a,  b,  and  o  .  However,  since 

0  <  <_  1,  it  follows  that  0  _<_  z  <  «,  and  hence  the  assumption  of  normality 

of  z  may  not  be  justified.  However,  the  least  square  technique  may  still  be 

2 

used  for  the  estimation  gf  a,  b,  and  a  .  Using  (26)  and  the  Taylor's  series 
estimation  procedure,  it  will  be  shown  how  one  can  obtain  estimates  of  Rx  and 
R„  from  the  estimates  of  a  and  b. 


All  these  quantities  can  be  numerically  computed.  The  normal  equations  (31) 
can  be  written  as 


a  all  +  b  a12  3  C1 
a  a2l  +  b  a22  =  c2 


Consider  the  matrix 


A  = 


11 

*21 


12 

S22 


Its  inverse  is 


,-l 


alla22"a12 


-a 


22 

12 


21 

'll 


Cl  early 


(3< 


(3f 


(3( 


and  the  estimates  a  and  b  of  a  and  b  are  Immediately  given  by  the  equations 


;  a22cl"a12c2 

a  =  __ 

a 1 la  22-a 12 

-a12c1+auc2 
b  =  - r  . 

alla22"a12 


(3 

(3! 


It  can  be  verified  that  a  and  b  are  unbiased  estimators  of  a  and  b  so  that 
E[a]  =  a  and  E[b]  =  b.  Knowing  a  and  b,  one  can  obtain  an  unbiased  estimate 

O 

of  a  (see  e.g.,  [4])  which  Is  given  by 

n 


aO  1  M  A  A  p 

0  *n X  (z1-aurbv1) 
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Further,  one  obtains  [5]: 


Var[a]  -  -i^Z  ;2 
aUa22_a12 

Var[b]  =  - g-  o2 

alla22'a12 

»  «  ai  p 

Cov[a,b]  = - — — j- 

alla22'a12 


or  more  succinctly,  the  covariance  matrix  of  (a,b)  is  given  by 
Var[a]  Cov[a,b]  \  *2  a22  -a12 

Cov[a,b]  Var[b]  J  aua22"ai2  “ai2  all 


(40) 

(41) 

(42) 


(43) 


b.  Estimates  of  R  and  R 

y 

It  is  clear  that  in  the  Taylor's  series  estimation  procedure,  the 
determination  of  ECP^]  and  VarfP^f]  requires  as  input  the  values  of 


Rx  •  E[Rx],  Ry  =  E[Ry] 
and  the  covariance  matrix 


(44) 


/  Var[Rx]  Cov[Rx,Ry]\ 

\  Cov[Rx,Ry]  Var[Ry]  / 


(45) 


The  intent  now  is  to  use  the  results  of  the  previous  least  square 
analysis  which  yielded  a  and  b,  the  estimates  of  a  and  b,  respectively,  as  well 
as 

j  Var[a]  Cov[a,b]  ' 

\  Cov[a,b]  Var[b] 
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to  obtain  (44)  and  (45).  This  can  be  accomplished  through  the  functional 
relationship  given  by  (26). 


Before  one  proceeds  any  further.  It  should  be  noted  that  the  notation  a 
is  used  to  denote  both  the  estimator  of  a  as  a  random  variable  and  the 
particular  numerical  estimate  of  a  obtained  from  the  least  square  analysis  and 
given  by  (37).  The  same  holds  true  for  b.  In  the  sequel,  the  meaning  of  the 
notation  a  or  b  should  be  clear  from  the  context. 

Since  a  and  b  are  not  known  exactly  but  are,  in  general,  random  variables 
it  follows  from  (26)  that  Rx  and  Ry  are  themselves  random  variables  given, 
respectively,  by 


r  =  -1  and  R  =  -! 

£  y  £ 


(46) 


Expand  now  Rx  =  a-1/2  as  a  Taylor's  series  about  the  point  a.  One  obtai 


ns 
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-  dRx 
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.  i  .v2 d  Rx 
-  ♦  jr  va-a)  — r* 

a=a  da 


a=a 


_1 

/a 


(a-a) 


(a-a)2  . 


(47) 


Recall  that  a  is  an  unbiased  estimator  of  a,  i.e.,  E[a]  =  a.  Taking 
expectations  on  both  sides  of  (47)  yields: 

\  ■  E<V  *4-^377  Et<a-a)2] 

=  4  +  —W  Var[a“]  •  (48) 

/a  8a3' ^ 

Since  a  is  an  unknown  quantity,  its  numerical  estimate  a  given  by  (37)  is  used 

a 

instead.  The  value  of  Var[a]  is  given  by  (40). 


yv 


Subtracting  (48)  from  (47)  one  obtains  the  following  up  to  and  including 


first  order  terms: 


Rx  ■  \  “  77*77  (a"a)  * 


Squaring  both  sides  of  (49)  and  taking  expectations  yields 


Var[R  ]  *  -i-  Var[a]  . 
x  4a  ^ 


Similarly  one  obtains 


*y  -  ~~  +  — It*  Var[b] 

y  rs  8bbU 


V*y  •  ^57?  (b"b) 


Var[R  ]  =  — j  Var[b]  . 
*  4b  ° 


Multiplying  (49)  by  (52)  and  taking  expectation  yields 


Cov[R  ,R  ]  * - ~TT7  Cov[a,b]  . 

x  y  4(ab)J/^ 


Expressions  (44)  and  (45)  are  thus  completely  specified.  The  following  is  a 


summary  of  the  results  obtained  when  a  and  b  are  replaced,  respectively,  by  their 


estimate  a  and  b 


R  =  —  +  — Var[a],  IT  =  —  +  — ir*  Var[b] 

x  /:  8a5/?  y  £ 

a  b 


Var[Rx]  Cov[Rx,Ry] 


Cov[Rx,Ry]  Var[Ry] 
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Cov[a ,b] 
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Another  Linear  Multiple  Regression  Scheme:  Method  2 


a 


This  section  digresses  on  a  second  linear  multiple  regression  scheme  to 
determine  estimates  of  the  parameters  Rx  and  Ry  of  the  Carleton  damage 


function 


\f  *  *XP  [- Cr?  *  *7)] 

Rx  Ry 


In  this  instance,  it  is  assuned  that  one  knows,  a  priori,  values  Rv  and  Ru 

*o  *o 

close  to  the  estimates  of  Rx  and  Ry,  respectively.  Specific  procedures  to 

obtain  such  initial  values  of  R  and  R  will  be  discussed  later. 

xo  yo 

The  expression  for  Pkf  is  expanded  as  a  Taylor's  series  about  the  point 
(Ry  .Ru  )  to  yield  up  to  its  first  order  terms 

An  n 


3P, 


kf 


kf 


kf 


r  r  +  ^x”^x  ^  aR 
»%  0  x 
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kf 


R,  ,Ry  *  ‘W  ^ 
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u  = 


2x 


Equation  (58)  becomes: 


z  =  au  +  bv  . 


RyRy0 


(60) 

(61) 


(62) 


a.  The  Least  Square  Technique 

Clearly  z  is  a  linear  function  of  u  and  v.  The  method  just  described 

approximates  what  would  have  resulted  in  a  nonlinear  regression  scheme  with  a 

linear  regression  scheme  using  the  Taylor's  series  expansion  and  retaining 

only  first  order  terms.  The  variables  u  and  v  are  obtained  from  the  known 

coordinates  (x,y)  at  which  is  measured  experimentally.  The  numerical 

values  of  z  at  a  particular  (u,v)  are  deduced  from  the  known  values  of  x,y  and 

Pk|r.  Under  these  conditions. one  is  again  dealing  with  a  linear  multiple 

regression  scheme,  and  the  values  of  a  and  b,  and  hence  R  and  R  ,  can  be 

*  y 

estimated  using  least  square  techniques  as  described  in  the  previous  section. 
The  output  results  of  the  least  square  approach  will  consist  of  the 

A  A, 

estimates  a  and  b  of  a  and  b*  respectively,  as  well  as  the  covariance  matrix 
Var[a]  Cov[a,b]  ' 

.  (63) 

Cov[a,b]  Var[b] 

From  (60)  it  follows  that 
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and,  of  course. 


(64) 
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It  is  possible  to  obtain  better  approximations  for  the  estimates  in  (64) 
and  (65)  by  retaining  higher  order  terms  in  the  Taylor's  series  expansion. 
This,  however,  would  considerably  add  to  the  complexity  of  the  problem.  A  more 


natural  approach  consists  in  using  the  initial  values  of  R„  and  Ru  to  obtain 

*o  yo 

estimates  Rx  and  R  using  the  linear  regression  scheme.  This  could  be 
1  “1 

followed  by  a  second  iteration  in  which  Ry  and  Rv  would  be  used  as  initial 

*1  *1 

values  to  obtain  improved  estimates  R„  and  Ru  using  the  same  linear 

*2  y2 

regression  scheme  and  so  forth,  assuming  that  the  procedure  will  converge  to 
some  value  of  R  and  R„. 

It  should  be  noted  that  one  of  the  advantages  of  the  present  multiple 
regression  scheme  over  Method  1  discussed  in  Section  IV-2  is  that  values  of 
pkf  =  0  need  not  be  discarded  from  the  data  in  the  least  square  analysis, 
provided  that  they  be  counted  as  additional  observation  points. 
t>.  Initialization  of  the  values  of  R.  and  R 


The  previous  method  on  the  linearization  of  a  non-linear  multiple 

regression  scheme  using  Taylor's  series  expansion  procedure  is  predicated  on 

the  assumption  that  initial  values  Rv  and  R„  for  the  two  weapon  radii,  which 

*o  *0 

are  close  to  the  estimated  values,  can  be  obtained.  Next,  three  possible 

avenues  are  discussed  that  could  be  used  to  obtain  such  initial  values  using 

experimental  data.  Such  data  is  usually  in  the  form  of  a  matrix  giving  values 

of  for  various  values  of  (x,y)  when  the  weapon  detonates  at  (0,0) » 

In  the  first  procedure,  the  values  of  R„  and  Rv  are  obtained  from  a 

*o  yo 

knowledge  of  the  mean  area  of  effectiveness  (MAE)  of  the  weapon  as  well  as  the 
impact  angle  I.  The  following  relation  is  used: 


MAE  =  /  /  Pkf  (x,y)  dx  dy  =  ^  R  . 


(66) 


The  value  of  the  double  integral  is  obtained  from  field  data  giving  Pk^(x,y) 

at  specific  values  of  x  and  y.  Assume  now  that  FL  is  the  weapon  radius  in 

* o 

the  direction  of  range  and  R  is  the  weapon  radius  in  the  direction  of 

yo 

deflection.  Let  I  be  the  weapon  impact  angle.  Then  the  ratio  R  /R 

xo  yo 

satisfies  the  following  experimental  relation  (see  e.g.,  [3]): 


1-0.8  cos  I  . 


(67) 


Using  (66)  and  (67)  values  of  Ry  and  R  can  be  obtained. 

o  ^o 

In  the  second  procedure  one  uses  the  following  relations  to  obtain  Ry 

*o 

and  R  .respectively,  (see  e.g.,  [6]) 

Jfo 


K  =  /  2x  P  f  (x  ,0)  dx 
o  0 


Ry  /  2y  Pkf  <°-y)  dy  * 


2 

Thus  R  is  obtained  by  multiplying  each  value  of  the  damage  function  on  the 
o 

x-axis  by  2x  and  integrating  the  result  over  all  possible  values  of  x. 

2 

Similarly,  Ry  is  obtained  by  multiplying  each  value  of  the  damage  function 
yo 

along  the  y-axis  by  2y  and  integrating  it  over  all  possible  values  of  y. 

To  use  the  third  procedure,  one  notes  that  in  the  expression  for  the 
Carleton  damage  function,  the  equiprobability  contour  line  corresponding  to  a 
probability  of  kill  of  P^f  =  e"*  =  .3678794,  has  for  equation 

2  2 

+  y  -  i 

2  +  ~7  ~  1  • 

R  R^ 

x  y 


Thus  the  two  weapon  radii  Rx  and  Ry  correspond  to  the  semi -axes  of  an  ellipse, 
and  these  values  can  be  determined  from  the  Pkf  matrix  obtained  from  the 
fragmentation  field  data. 


SECTION  V 


CONCLUSIONS 

It  has  been  customary  to  evaluate  weapon  effectiveness  by  the  use  of  a 
single  index  number  such  as  the  probability  of  kill  P^.  In  the  present  report 
procedures  have  been  developed  for  estimating  the  error  in  the  P^  value 
knowing  the  error  in  the  input  parameters  used  to  compute  Pk. 

For  fragment  sensitive  targets,  the  Carleton  damage  function  appears  to 
be  a  reasonable  model  to  compute  P^  in  the  absence  of  blast  and  aiming 
error.  This  function  usually  contains  the  two  parameters  Rx  and  Ry  which 
define  the  weapon  radii.  The  report  has  addressed  itself  to  the  specific 
problem  of  utilizing  the  inherent  errors  in  Rx  and  Ry  to  compute  the  error  in 
pk- 

Methodologies  have  been  developed  to  identify  the  causes  of  error  in  the 
parameters  when  field  data  are  used  to  estimate  Rx  and  Ry.  These 
methodologies  are  based  on  well  established  statistical  techniques  such  as 
linear  multiple  regression,  and  their  use  has  provided  estimates  for  R  and  R 

X  j 

as  well  as  estimates  for  Var[Rx],  Var[Ry], and  Cov[Rx,Ry].  In  turn,  these 
estimates  can  be  used  to  compute  E[Pk3  and  Var[Pk].  Knowing  the  first  two 
moments  of  ,  Chebyshev's  inequality  provides  one  with  a  means  for  setting  up 
confidence  intervals  for  Pl.  values. 
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